Antiferromagnetic vs ferromagnetic interactions and spin-glass-like behavior in 

ruthenates 
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We have made a series of gradient-corrected relativistic full-potential density-functional calcula- 
tion for Ca-substituted and hole-doped SrRu03 in para, ferro, and A-, C-, and G-type antiferro- 
magnetic states. Magnetic phase-diagram data for Sri_^Ca^Ru03 at K are presented. Neutron 
diffraction measurement combined with total energy calculations show that spin-glass behavior with 
short-range antiferromagnetic interactions rules in CaRuC>3. The substitution of Sr by Ca in SrRu03 
decreases the ferromagnetic interaction and enhances the G-type antiferromagnetic interaction; the 
G-AF state is found to stabilize around x — 0.75 consistent with experimental observations. In- 
clusion of spin-orbit coupling is found to be important in order to arrive at the correct magnetic 
ground state in ruthenates. 

PACS numbers: 75., 81.05.Je, 75.10.Nr, 71.20.-b 
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Ever since unconventional superconductivity was ob- 
served in Sr2RuC>4 Q], ruthenates have attracted much 
interest. Substitution || of the smaller Ca for Sr in 
Sr 2 Ru04 leads to an antiferromagnetic (AF) Mott in- 
sulator with a staggered moment of S = 1. Coexistence 
of ferro- and antiferromagnetic fluctuations in Ca2Ru04 
and competition between p- and d-wave superconductiv- 
ity in Sr2RuC>4 have also been reported M. Recently 
coexistence of magnetism and superconductivity is dis- 
covered J! in ruthenium-based layered cuprates such as 
i?RuSr 2 Cu 2 8 (R = Eu, Gd). Lee et al. fj reported 
the existence of a pseudogap (PG) in BaRu03 remi- 
niscent of the PG in high-T c superconductors. As in 
/-electron-based intermetallics, non-Fermi-liquid (NFL) 
behavior has recently been observed in La4RugOi9 ||. 
Critical magnetic fluctuations associated with metamag- 
netism is reported for Sr3Ru207 Q] and a metallic AF 
phase with temperature induced insulator-to-metal tran- 
sition is found for Ca3Ru207 ||. 

Poorly metallic NFL-behaving SrRu03 is the only 
known ferromagnetic (F; Tq = 160 K) 4d transition-metal 
oxide || |l0|. CaRu03 is also metallic, but experimen- 
tal and theoretical studies conclude contradictory regard- 
ing the nature of the magnetic ground state (e.g., AF 
H |l2|| nearly F |r|, exchange enhanced paramag- 
netic (P) Q, Curie- Weiss P verge of F instability 
Jill , spin glass (SG) etc.), and this controversy is 
not settled yet. Some of these reports indicate a lack 
of long-range magnetic order whereas others suggest ev- 
idence for an AF ground state with a Neel temperature 
(Tjv) of -110K Q. The striking difference in the mag- 
netic properties of SrRuOs and CaRuOs makes magnetic 
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phase diagram studies on Sri-^Ca^RuOs quite interest- 
ing. 

CaRuC>3 and SrRuOs are isostructural and isoelec- 
tronic, differing structurally only in the degree of the 
small orthorhombic distortion (lattice constants within 
2%) . Hence, insight into the magnetic properties of these 
compounds is expected to increase the general knowledge 
on magnetic phenomena in perovskite oxides, and ruthen- 
ates in particular. The density-functional calculations 
]l6| , |l8| for CaRu03 which conclude with "verge of F in- 
stability" have the weakness that AF interactions were 
not considered. In this letter, we present the intriguing 
result that SG behavior with short-range AF interaction 
occurs in metallic CaRuO^. 

The present full-potential linear-muffin-tin orbital [Tg| ] 
(FPLMTO) calculations are all-electron, and no shape 
approximation to the charge density or potential has been 
used. The basis functions, charge density and potential 
were expanded in spherical harmonic series inside the 
muffin-tin spheres and in Fourier series in the interstitial 
regions. The calculations are based on the generalized- 
gradient-corrected (GGA) density-functional theory as 
proposed by Perdew et al. pp| . Spin-orbit (SO) terms 
arc included directly in the Hamiltonian matrix elements 
for the part inside the muffin-tin spheres. The basis set 
contained semi-core 4p and valence 5s, bp, and 4d states 
for Sr, 4s, 4p, and 3d states for Ca, 5s, 5p, and Ad for Ru, 
2s, 2p, and 3c? states for O, and 4s, 4p, and 3c? states for 
K. All orbitals were contained in the same energy panel. 
The self consistency was obtained with 284 k points in 
the irreducible part of the first Brillouin zone for the or- 
thorhombic structure and the same density of k points 
were used for the cubic structure as well as for the super- 
cells. Structural parameters for Sri^Ca^RuOs are taken 
from the X-ray diffraction data reported by Kobayashi et 
al. pl[ . Details about the experimental reexamination 
of CaRuOs will be presented in the full account of our 



2 



TABLE I: Total energy (relative to the lowest energy state; in 
meV/f.u.) for CaRu0 3 in P, F, and A-, C-, and G-AF states 
with orthorhombic and undistorted cubic (designated cubic) 
structures. 
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C-AF 


G-AF 
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findings 

There are three possible magnetic arrangements ac- 
cording to the interplane and intraplane couplings in per- 
ovskite oxides, (i) With interplane AF coupling and in- 
traplane F coupling the A-AF structure arises, (ii) The 
opposite arrangement with interplane F coupling and in- 
traplane AF coupling is called C-AF structure, (iii) If 
both the inter- and intraplane couplings are AF the G- 
AF structure arises. In the G-type AF lattice, each Ru 
atom will be surrounded by six Ru neighbors whose mo- 
ments are antiparallel to that of the chosen reference 
atom. The details about the calculations on magnetic 
and excited state properties can be found elsewhere p3| . 
We have placed the magnetic moment direction along 
[001] in all calculations to comply with the out-of-plane 
easy axis observed for Sri_ x Ca K Ru03 Q. The energy 
differences between the various magnetic configurations 
are very small and extreme computational accuracy is 
therefore needed. The theoretical method used in the 
present study is capable of reproducing energy differences 
of the order of [ieV [^5| . 

The calculated total energy for CaRuC>3 from the GGA 
calculation with (GGA+SO) and without SO coupling- 
are given in Table [| relative to the corresponding ground 
state. From this table it is clear that the G-AF phase 
is the ground state for CaRu03. The total energy dif- 
ferences between the different AF phases are as already 
mentioned very small. If we disregard the orthorhom- 
bic distortion, the F phase is seen to have the lowest 
energy indicating that strong magneto-elastic couplings 
are present in CaRuOs. We also see large energy gains 
when we include the orthorhombic distortion in the cal- 
culation. The F phase is found to be the ground state 
with 7 meV/f.u. lower energy than the G-AF phase when 
SO interactions are neglected. Hence, SO coupling ap- 
pears to play an important role in deciding the magnetic 
properties of ruthenates. 

The calculated magnetic moment at the Ru site in 
SrRuOs is 0.92/ie/atom is in good agreement with 
0.85 /is/atom obtained from magnetization measure- 
ments Q. Goodenough |2(| has argued that, in general, 
a spontaneous change between F and AF states are ex- 
pected when hf ~ P^af, the stable magnetic phase being 
the one with the highest atomic moment. Our findings 
are consistent with this viewpoint in that the larger mag- 
netic moments are found for the stable magnetic configu- 
ration of Sri-zCa^RuOs (viz., for x < 0.724 the F phase 
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FIG. 1: Magnetic phase diagram data for Sri-^Ca^RuOa. 
Energies are given with respect to the total energy for the F 
phase. 



has highest magnetic moment and for x > 0.724 G-AF 
has the highest moment). The hyperfine field at the Ru 
site in G-AF state CaRu03 is 263kOe//is in agreement 
with the experimental value of 222 ± 50 kOe//is obtained 
from Knight shift measurements J27J . 

The experimentally found spontaneous magnetic mo- 
ment |l3| ] and Weiss temperature p8| ] for Sri_ x Ca^RuOg 
are given in Fig. |l| along with the calculated total energy 
difference relative to that of the F phase for the various 
magnetic configurations. 

An increase in x decreases the F interaction and in- 
creases the G-AF interaction. Because of the smaller 
ionic radius of Ca compared with Sr (0.99 vs 1.18 A), 
tilting of the octahedra occur when one substitute Ca for 
Sr in SrRuOs which in-turn affect the magnetic proper- 
ties. Fig. [j] shows that the F to G-AF transition takes 
place around x — 0.75 both experimentally and theoret- 
ically. This confirms the presence of AF interactions in 
CaRuOz. One of the interesting aspects of Fig. [I] is that 
at x = 0.5, there is large gain in the total energy of F 
phase compared with P, A- and G-AF phases. DOS anal- 
yses show that the F phase of Sro.sCao.sRuOs is nearly 
half-metallic with a gap of 0.73 eV in the majority spin 
channel which gives extra contribution to the stability. 
Moreover, Ep is at peaks in the DOS curves of the P, 
A- and G-AF phases (which further contribute to the 
enhanced energy difference relative to the F state). 

The lower total energy for the G-AF state for x > 
0.724 compared with the other magnetic configurations 
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FIG. 2: Total DOS for CaRuOa in F and G-AF states and 
hypothetical KRu03 in G-AF state. Fermi level is set to zero. 



can be understood as follows. The spin-projected total 
DOS curve for CaRu03 in the F and G-AF phases are 
shown in the lower and middle panel of Fig. ^. In the 
F case Ep is at a peak-like feature in the DOS curve. 
This infers an instability to the lattice. However, for the 
G-AF phase a deep valley-like (PG-like) feature appears 
in the DOS curve. This will give an extra contribution 
to lattice stability p3] and hence CaRuOa is expected to 
stabilize in the G-AF phase. The calculated excited-state 
properties such as reflectivity, O K-edge spectra and XPS 
spectra for CaRu03 in the G-AF phase are found to be in 
good agreement with experimental spectra (see Fig. ||), 
indicating that the ground state is correctly assigned. The 
creation of a PG-like feature on introduction of G-AF 
ordering gives an extra contribution to the band energy 
term of the total energy (66meV/f.u.) which stabilizes 
the G-AF phase over the F phase. However, our neu- 
tron powder diffraction diagrams at 298 and 8 K show 
no sign of extra Bragg reflections at 8K, which should 
have been the proof of long-range AF ordering. (A full 
account of the neutron diffraction findings will be given 
in the forthcoming article Hence, the magnetic ex- 

change interactions must have only short-range influence. 
Closer examination of the DOS curve for the G-AF phase 
shows that Ep is located on a shoulder of a peak-like fea- 
ture (an unfavorable condition for stability) with a local 
valley (LV, see Fig.||) situated just 40meV above Ep. 
This may explain why there is no established long-range 
magnetic ordering in this material. Recent magnetiza- 
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FIG. 3: Calculated valence band XPS, O K-edge XANES, 
and optical reflectivity spectrum for G-AF state CaRuC>3. 
Experimental XPS, XANES and reflectivity spectra are taken 
from Refs. §(], |||. 



tion measurements show distinctions between in the 
zero-field-cooled and field-cooled characteristics that may 
indicate SG behavior. So, we conclude that CaRuOs ex- 
hibits short-range AF ordering which manifests itself in 
SG behavior at low temperature. 

Rigid-band-filling analysis show that an addition of 
0.136 electrons/f.u. will bring Ep to LV with the ex- 
pected effect of long-range G-AF ordering. In fact, 
long-range magnetic ordering has been observed for 
CaRui-zRh^Oa with x > 0.015 where Rh acts as an 
electron donor for CaRuOa |3|] . Rigid-band- filling anal- 
ysis also show that removal of exactly one electron from 
CaRuOs will bring Ep to a PG. In order to test this 
possibility we have made additional calculations for hy- 
pothetical KRUO3 with the structural parameters for 
CaRuOs. Consistent with the rigid-band- filling findings 
Ep is found at PG for the G-AF phase of hypotheti- 
cal KRUO3 (see upper panel of Fig. g). When Ep is at 
PG, all bonding states would be filled and all antibond- 
ing states empty, implying extra contributions to stabil- 
ity. Total energy calculations show that the G-AF phase 
is 21.7meV/f.u. lower in energy than the F phase for 
hypothetical KRUO3, viz. a larger difference than that 
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for CaRu03 indicating that one really should expect G- 
AF long-range ordering in KRUO3 at low temperature. 
Further the G-AF phase is 38.6 meV/f.u. lower in en- 
ergy than the P phase for hypothetical KRUO3 indicat- 
ing that Tjy should be reasonably high. It is worth to 
recall that according to Goodenough |3j| the sign of the 
transfer integral for 180° cation-anion-cation superex- 
change interactions between octahedral-site cations are 
predicted to give AF for Ru 5+ -0-Ru 5+ . Ru 5+ ions with 
enhanced AF interaction has recently been established 
for CaRuo.95Cuo.05O3 Q. Moreover, the change in the 
oxidation state of Si'3Ru2 07 from Ru 4+ to Ru 5+ by fluo- 
rine addition stabilizes G-AF ordering ]35| . Oxygen non- 
stochiometry and hole doping can accordingly bring Ef 
toward PG and stabilize the G-AF ordering. This may 
explain why long-range AF ordering has been reported 
for some experimentally studied CaRuOs samples Jll[] 
and Na doped CaRu0 3 §§]. 

In the SG state the moments are arranged in a certain 
equilibrium oriented pattern, but without long-range or- 
der. Some characteristics of SG materials are: (i) Ex- 
istence of local moment as we have found with spin- 
polarized calculations, (ii) No magnetic Bragg scatter- 
ing at low temperature as our neutron diffraction study 
show, (iii) A history-dependent magnetic response as 
the recent magnetization measurements |l7j have shown, 
(iv) Both long-range and short-range terms in the Hamil- 
tonian p7j |. In accordance with our calculations we be- 
lieve that short-range nearest-neighbor AF interactions 
are dominant in CaRu03. 

Two important ingredients necessary to produce SG 
behavior are frustration and partial randomness of the 
interaction between the magnetic moments [ f38[ . Owing 
to the small total energy differences between the AF and 
F phases of CaRuOs (see Table| and Fig. [|), there will 
be a competition among the different AF and F inter- 
actions, in the sense that no single configuration of the 
the moments is uniquely favored by all interactions {viz. 
frustration). In other words, the AF moments can ar- 
range themselves randomly in small domain-like regions 
with minimal loss in energy. The Ru ions which are re- 
sponsible for the magnetic properties of CaRuOs are in 
a strictly periodic order. Hence, structural disorder can 
not be responsible for the SG behavior. The SG state 
of CaRuOs is characterized by a predominant AF situa- 
tion. AF systems with SG-like transition have been ob- 
served in ruthenates with pyrochlore-like structures such 
as -R2RU2O7 for which recent experimental results indi- 
cate that the atomic arrangement does not participate in 
the P-to-SG transition which is solely associated with the 
Ru moments f^] . The small energy differences between 
the different AF states bring disordering of the moments 
and hence SG behavior to CaRuOs at low temperature. 
In fact, SG behavior has been found for CaRuOs doped 

or Rh [B2| . The calculated 



with small amounts of Sn 1 40 



differ significantly between the SG and G-AF phases. 

We conclude that CaRuOs exhibits short-range AF in- 
teraction with SG behavior. The hitherto hypothetical 
KRUO3 is expected to exhibit long-range G-AF ordering. 
We have demonstrated that the relative strengths of the 
F and AF exchange interactions can be varied by varying 
the Ca content of Sri_ x Ca 2: Ru03. 

PR is grateful to the Research Council of Norway for 
financial support. Part of these calculations were carried 
out on the Norwegian supercomputer facilities. 



excited state properties for CaRuOs in the G-AF state 
are in good agreement with the experimental measure- 
ments indicating that the electronic structure does not 
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